Quantum correlations in topological quantum phase transitions 
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We study the quantum correlations in a 2D system that possesses a topological quantum phase 
transition. The quantumness of two-body correlations is measured by quantum discord. We calcu- 
late both the correlation of two local spins and that of an arbitrary spin with the rest of the lattice. 
It is notable that local spins are classically correlated, while the quantum correlation is hidden in the 
global lattice. This is different from other systems which are not topologically orderd. Moreover, the 
mutual information and global quantum discord show critical behavior in the topological quantum 
phase transition. 
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I. INTRODUCTION 

Topological phase is a new kind of order that cannot be 
described by the symmetry-breaking theory [1]. A typi- 
cal example is the quantum Hall system, which exhibits 
a lot of amazing properties, such as topological degener- 
acy and fractional statistical behaviors . Especially, the 
property of topological protection may lead a new way 
for quantum computation [2]. 

Different from the quantum Hall system, Kitaev toric 
code model is an exactly solvable spin lattice model that 
is topologically ordered [3]. The system is immune to 
small perturbations. The breaking down of the topo- 
logical phase happens through a quantum phase tran- 
sition [4]. A lot has been studied about the topological 
quantum phase transition, especially about the toric code 
model in the present of a magnetic field [5, 6]. 

Concepts of quantum information are borrowed to the 
study of quantum phase transition, like entanglement and 
fidelity [7, 8]. Here, we are interested in the correlations 
in topological phase , because the magic power of quan- 
tum computation roots from the strange non-classical 
correlations. 

Entanglement is the most important non-classical cor- 
relation in quantum information processing, such as 
quantum teleportation [9]. However, some separable 
states also have properties that are not achievable by 
classical methods [10]. Recent results suggest that these 
correlations may also take effect in quantum computa- 
tion. The classification of non-classical correlations and 
their effects still remains an open problem [11]. 

Quantum discord is a measurement for the "quantum- 
ness" of a pairwise correlation [10]. It is based on the 
fact that the mutual information has two equivalent def- 
initions in the classical world, while their quantum gen- 
erations are not equivalent. Quantum discord is defined 
as the minimum of their difference and measures how 
"quantum" the correlation is. Besides entangled states, 
some separable states also have non-zero quantum dis- 
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cord, which means they are non-classical. Recent studies 
suggest that quantum discord but not entanglement may 
be responsible for mixed-state computation [12, 13]. 

The quantum discord in quantum phase transition 
were studied in ID systems, such as XXZ-ch&in and 
some other ^-symmetric ID spin models [14, 15]. Dur- 
ing the phase transition, quantum discord shows criti- 
cal behavior at the phase transition point. The study 
of quantum discord in thermal Heisenberg system also 
shows some different behavior from entanglement [16]. 

In this paper, we study the correlations in Castelnovo- 
Chamon model [17], which is a 2D system that shows a 
quantum phase transition from a topologically ordered 
phase to a magnetized one. It is a deformation of toric 
code model and possesses higher symmetry than other 
ID models mentioned above. Both local spin-spin corre- 
lation and that between a spin and the rest of the whole 
lattice are calculated. It is notable that in such a topo- 
logically ordered system, quantum discord of local spins 
is always zero in both phases, which means the local cor- 
relations are totally classical. While the correlation be- 
tween a local spin and the rest of the lattice behaves 
more like a pairwise entangled pure state, and the quan- 
tum discord signals critical point. This is different from 
previous studies in other models [14, 15, 16]. Our results 
shows that in topologically ordered system, the quantum 
correlation is hidden in the lattice globally by the high 
symmetry of the system. 

Besides, we calculated the mutual information of the 
correlations in the system. It was pointed out that in 
topologically ordered systems, which have no local order 
parameter, the topological quantum phase transition can 
be signaled by local properties like the reduced fidelity of 
two spins, and it is even more sensitive than the global 
fidelity [18, 19]. Here, we calculate the mutual informa- 
tion, both of the global correlation and that of two local 
spins. We see that the mutual information could also 
characterize the critical behavior. 

The paper is organized as follows. In Sec. II we briefly 
review the concept of quantum discord and the basic 
properties. In Sec. Ill, we introduce the Castelnovo- 
Chamon model. We give the ground state and explain 
how it can be mapped to the classical Ising model. In 
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Sec. IV, we calculate the quantum discord of local spin- 
spin correlations, and also that of a local spin with the 
rest of the whole lattice. Conclusion is drawn in Sec. V. 



II. QUANTUM DISCORD 

Quantum discord can be used as a measure of the quan- 
tumness of a pairwise correlation [10]. In this part, we 
briefly introduce this concept. 

Information is amount of uncertainty that can be elim- 
inated after we get a measurement result. Mutual infor- 
mation X(A : B) describe the information about A we 
gain after the measurement of B, or rather, the infor- 
mation that A and B have in common [20]. In classical 
world, there are two equivalent definitions of mutual in- 
formation, 

1(A:B) = H(A) + H(B) - H(A, B), 
J{A : B) = H(A)-H(A\B), 

with the Shannon entropy H(-) = —J2iPi^°E2Pi- 
H(A\B) is the conditional entropy, which is a measure 
of how uncertain we are about A, averagely, when B 
is known. The Bayes' law tells us that p(ai,bj) = 
p(ai\bj)p(bj) = p(bj\di)p(ai) , where p(a,i\bj) is the con- 
ditional probability that describes the probability to get 
o, when we know the value of B is bj . That guarantees 
the equivalence of the two definitions above. 

However, things are different when we generalize the 
concepts to the quantum world. We can get the quan- 
tum version of 1{A : B) easily by replacing H with von 
Neumann entropy S(p) = — tr(plogp). While the con- 
cept of conditional entropy in fact implicitly calls for a 
measurement of B. To get the corresponding J (A : B), 
we have to choose a set of basis {Ilf } to measure sys- 
tem B. The state of the system after measurement is 
Pi = tlfpAstl-flPi, where Pi = tr(ftf PAB^f )■ With 
the knowledge we gain after the measurement, we get 



J(A\{flf}) = S( PA )-S( PAB \{uf}) 
= S(p A ) - ^2piS(pi). 



(1) 



The value of J{A\{Uf}) depends on the choice of {ftf }, 
i.e., how we measure B, and therefore may not be equal 
to J (A : B) any more. Quantum discord is defined as 
the minimum of their difference, 



D{p) = mm l(A:B)-J(A\{flf}) 



(2) 



^/(^{nf }) describes the amount of information of A 
achievable by projective measurements on B. It can be 
proved that D(p) > 0. From the derivation above, we see 
that the quantum discord of classical correlations should 
be zero. We can use quantum discord as a measure of 
the "quantumness" of a two-body correlation. 
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Figure 1: (Color online). Demonstration of the model. Spins 
lie on the edges (like The nearest spins become 

next-nearest in the dual lattice (the dashed line). The red 
plaquette and the blue cross represent the plaquette and star 
operators. The red line across the lattices represents a prod- 
uct of of along the non-trivial loop on a torus. The system 
is invariant under this transformation. 



The quantum discord of states that contain entangle- 
ment is obviously non-zero. It should be emphasized 
again that not all separable states are classical under the 
definition of quantum discord. A simple example is p = 
(|00)(00| + I + +)(++ I) /2, where |+) = (|0) + |1))/V2. 
The quantum discord is non-zero. The information inside 
the state cannot be fully extract just by local projective 
measurements. From the description above we can see 
that the quantum discord of pab is zero, if and only if 
Pab has the form of pab ~ J^PiPt ® nf '. 



III. CASTELNOVO-CHAMON MODEL 

Castelnovo and Chamon proposed a model that shows 
topological quantum phase transition [17]. It is a defor- 
mation of Kitaev toric code model. The Hamiltonian is 



H — — \q 2^ Bp — Ai 2^ As 




and B 



where Ao,i > 0, A s = J7 

the star and plaquette operators in toric code model re- 
spectively. (3 is a coupling constant. The star operator 
A s acts on the four spins around the vertex s, while the 
plaquette operator B p acts on the four spins on the edges 
of the plaquette q, as shown in Fig. 1. We consider the 
problem under the torus boundary condition. 

The ground state can be written down analytically. We 
give the state in the topological sector that contains the 
fully magnetized state |0) = | |TT ■ ■ ■ T) as 



|GS(/3)) = z(py 
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S|0>, (4) 



with Z(fi) = E ffG G ex P[' 9 Ei a 'i (S 1 )]- G is the Abelian 
group generated by the star operators {^4 S }- So g\0) con- 
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tains separable spins taking the form like |011 1 10 . . . 0) , 
and we can denote each g\0) by a corresponding binary 
number as \x) (as what we do in the following). crf(g) 
is the value of spin at site i in state g\0). The sum in 
the exponential term in fact counts the total magnetic 
polarization of g\0). 

It may be not obvious to get Eq. (4) directly. We can 
just put it back into Eq. (3) and it can be easily checked. 
When (3 = the model reduces to the toric code model. 
When — > oo the ground state becomes the fully magne- 
tized reference state |0). At j3 c = (l/2)ln(\/2 + 1) there 
is a second-order topological quantum phase transition, 
according to the study of topological entropy [17] and 
fidelity [18, 19] in this model. 

Furthermore, the value of of (g) (note that there is no 
hat, it is just an integer number relating to g) of the 
ith spin is actually determined by whether the two ends 
(s and s' in Fig. 1) of the zth edge are acted by A s ( s /) 
or not. As Aj. = 1, elements of G can be represented 
as a configuration of {0 S }, where 6 S = +1 means A s 
acts on vertex s, while 6 S = — 1 means not. So we get 
°f (fl 1 ) = 9s9s'- The normalizer in the ground state Eq. 
(4) is 



z{p) = cx p 



p Y 9 > 9 > 



(5) 



which is just the canonical partition function of 2D clas- 
sical Ising model without external field, with the Hamil- 
tonian ifi sing = - Yl( ss >) Ms'- 

As an example, we can calculate the correlation func- 
tion as 




Figure 2: Mutual information of two nearest spins of of// (the 
solid line) and crfcr*, (the dashed line) as shown in Fig. 1. 
Critical change happens at /3 C = (l/2)ln(A/2 + 1), which is 
in accord with previous studies. The quantum discord of the 
two spins are always zero. 



A. Local spin-spin correlation 

Firstly, let us look at the correlation of two local spins. 
We need to get the reduced density matrix of two spins. 
The set where p, v take value 0,1,2,3 and 

<7° = 1, contains 16 matrices and they form a complete 
orthonormal basis for 4x4 Hermitian matrices under the 
Hilbert-Schmidt inner product (A,B) K . S = tt(A*B) [20]. 
Conveniently, the reduced density matrix can be wirtten 
as the expansion of the basis set {^fo-j} [18, 22], 



Pij 



3 



(7) 



(GS|<t?|GS) = ]>>*(<?) exp 
= Y 9 s e s > exp 

= (#0,0#0,l)lsing- 



p Y 

(as') 



(6) 



Thus, we can see that the model can be mapped to the 
Ising model, which is exactly solvable [21]. 



where (<rf ) = Tr(p<3gd'| i <7j) = tr(p ij (rf <rj) is the inner 
product of pij and <rf <rj. 

Furthermore, most terms above can be eliminated be- 
cause of the symmetry of the system. Draw a closed 
loop through the torus arbitrarily (as the red line shown 
in Fig. 1), and define a corresponding transformation 
V = Ylu ne of . The Hamiltonian Eq. (3) is invariant 
under the transformation V. Also, should commute 
with any V. Only the terms l,of,<7| and df<r| could 
exist, so we get 

hi = + + + (8) 



IV. CORRELATIONS IN THE LATTICE 

In this section, we discuss the correlations in the lat- 
tices. Both the local spin-spin correlation and the global 
correlation between a single spin with the rest of the lat- 
tice are considered. We find that the local correlations 
are classical, and the quantum correlation emerges only 
when considering the whole lattice. In both cases, the 
mutual information signals the critical behavior. 



The density matrix is diagonal. It can be written in the 
form of pij ~ Pnpn <8> n„ . According to what we have 
seen in Sec. II, the quantum discord of pij is zero. That 
means the correlations between any two local spins are 
always classical. This is quite different from other studies 
of quantum discord in the phase transition of ID systems 
that are not topologically ordered [14, 15, 16], where the 
quantum discord of local spins shows different behavior 
in different phase areas and exhibits critical behavior. 
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This roots from the high symmetry of the topologically 
ordered system. This 2D system exhibits higher symme- 
try than other ID ^-symmetric models [14, 15]. The 
system is conserved under the transformation of V along 
any closed loop, which eliminates all non-diagonal terms. 
So the quantum discord is zero in the topological phase 
area. It was stated in Ref. [2] that in a topologically 
ordered system, all observable properties should be in- 
varint under smooth deformations (diffeomorphisms) of 
the space-time manifold, which means the only local op- 
erator that has nonvanishing correlation functions is the 
identity. For example, in toric code model, only iden- 
tity exists in the expansion Eq. (7) and /3y ~ 1 ® 1, 
which means local spins are even uncorrelated. While 
the other phase area where (3 — > oo is a fully magnetized 
phase, which obviously only contains classical informa- 
tion of probability . This is why local correlations are 
classical in both phases. 

Besides, it was stated that topological quantum phase 
transition cannot be described by the symmetry-breaking 
of a local order parameter and involves a global rear- 
rangement of non-local correlations [1]. However, re- 
cent researches indicated that some concepts in quan- 
tum information theory, which describe local properties 
although, still signal the singularity in topological quan- 
tum phase transition [5, 17, 18]. The reduced fidelity 
and local magnetization were studied in the Castelnovo- 
Chamon model and the Kitaev toric code in a magnetic 
field and they exhibit critical behavior of the topologi- 
cal quantum phase transition. We calculate the mutual 
information X of two nearest spins, which is also a lo- 
cal property (Fig. 2). The correlations in pij can be 
evaluated with the help of the mapping to Ising model, 
as mentioned in Sec. III. We can see that the mutual 
information of both nearest and next-nearest spins (in 
the dual lattice) exhibits critical behavior. But the next- 
nearest mutual information is much less sensitive. 



B. Global correlation in the lattice 

As we have seen in the last part, the correlations be- 
tween local spins are completely classical in both phases. 
In this part, we calculate the correlation between a local 
spin and the rest of the whole lattice. As the increasing 
of /3, the system turns to the magnetic phase, and the 
correlation between a local spin and the lattice becomes 
more and more "classical". 

To calculate the correlation of an arbitrary spin de- 
noted by k with the lattice, we treat the rest of the lat- 
tice as a whole system. We can always rewrite the ground 
state as 

|GS(/3)) = ^(bKlOH + JWIl)* 



a\X)\0) k + b\Y)\l) k . 



(9) 



where \X) = ^2 x a x \x) and \Y) = J2 y b y\y)- Tne x iV 
in the basis vectors are the binary number representa- 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



Figure 3: Global correlations between a local spin and the rest 
of whole lattice. Here, the quantum discord (the solid line) is 
equal to the entropy of entanglement, and just one-half of the 
mutual information (the dashed line) of the pairwise system. 
Both the quantum discord and the mutual information show 
critical change at the phase transition point. 



tion of g\0) excluding the kth spin, as mentioned in Sec. 
III. Notice that g\Q) and g'\0){g ^ g') have at least four 
different spins. So we are sure that \X) and \Y) have 
no term in common, and (A|y) = 0. Therefore, we can 
treat |GS(/3)) as a simple 2x2 entangled state. In this 
case, the quantum discord is equal to the entanglement 
of entropy [23, 24], 



D( PAB ) = 1(A : B)/2 = S(A) = S(B). 



(10) 



We calculate it in detail. The coefficients a x ,b y are 
superposition coefficients in Eq. (4) correspoindingly. So 
the value of a 2 (b 2 ) is just the Ising partition function with 
a constraint that cr| = 1(— 1), in another word, 9 r 6 r > = 
1(— 1), where r and r' are the nearest vertices of spin k. 



12 exp 



{e e },e r e r , 

* = E 

{s B },e r e r , 



exp 



J2 9 ° d * 

(9s8 s ,) 



/Z{(3). (11) 



Notice that a 2 — b 2 is just the nearest correlation func- 
tion (0o,o0o,i)> 



70,0C0,1/ — 



e r e r ,=i o r e r ,=-i 
= a 2 -b 2 . 



(12) 



Together with a 2 + b 2 = 1, we can get the value of a, b. 

Now we calculate the quantum discord of the ground 
state Eq. (9). Instead of doing all the possible projective 
measurement to the spin, equivalently, we implement all 
possible local unitary operations on the spin and then 
measure it by {|0)(0|, |1)(1|}. 
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n tfp QS uii = p ®\o){o\ 

( a 2 cos 2 § ±ab sin 6e i(p \ 
~ 6 2 sin 2 § 

niiTtjSagEm^piiaiixii 

/ a 2 sin 2 f -ia6sin6»e^\ 
~ y-ldbsmee-** 6 2 cos 2 f ^® Ul ' 

where 

/ cosf sinfe-^A 
^-(.sinl^ cosf J' (13) 

The unnormalized post-measurement density matrices 
po and pi both have only one non-zero eigenvalue re- 
spectively, i.e., A fc = [a 2 + b 2 + (-l) k (a 2 - 6 2 )cos6l] /2, 
where k = 0, 1. That means the conditional information 
about the lattice after the measurement of a local spin is 
zero, and J is always equal to the entanglement of en- 
tropy S, no matter what measurement we impose on the 
local spin. So the quantum discord is equal to S, 

D(pgs) = J=1/2 = S 

= -a 2 log a 2 -6 2 log& 2 , (14) 

where 

a 2 = (1 + <0o,o0o.i» A 

b 2 = (1 - <0 o ,o0o,i» /2. (15) 

The quantum discord and mutual information of the 
global correlation is shown in Fig. 3. Comparing with 
that of local spins correlation, the quantum discord is not 
zero, which means the quantum correlation exists in the 
lattice globally. It also signals the critical point in the 
phase transition, just like the mutual information. As 
the increase of /?, the quantum discord decreases to zero, 
which means the global quantum correlation disappears 
gradually . 

In summary, the quantum correlation hides in the 
global lattice. We can only get classical correlations be- 
tween local spins. All these results of correlations seem 



to suggest that the ground state of the topologically or- 
dered system behaves like a generalized GHZ state. The 
quantum information is encoded in the lattice globally 
and so it can be protected better than in other systems. 



V. CONCLUSION 

In this paper, we studied the correlations in 
Castelnovo-Chamon model. Both local and global corre- 
lations were studied. The correlations were measured by 
quantum discord. As we have seen, local spins are clas- 
sically correlated although the Hamiltonian is so com- 
plicated. While the quantum correlation is hidden in 
the lattice globally. This is quite peculiar comparing 
with previous studies. We analyzed that these distinctive 
characters result from the high symmetry of the 2D topo- 
logically ordered system. The spins along any loop on the 
torus are ^-symmetric. This strict constraint clears the 
quantum correlations between local spins. Only global 
quantum correlation exists, just like a generalized GHZ 
state. We believe that this is a generic property in topo- 
logical quantum phase transition because of the partic- 
ular symmetry of topologically order systems, as men- 
tioned previously. 

Moreover, we calculate the mutual information of two 
nearest spins, which signals critical behavior of the topo- 
logical quantum phase transition. Similar to previous 
study of fidelity, mutual information also works as a lo- 
cal probe of the topologically ordered phase, although 
topological order cannot be described by the symmetry- 
breaking of local order parameter. More study is required 
for the correlations in other more realistic systems. 
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